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Abstract 

Neutrinoless double-/? decay is of fundamental importance for the determining neutrino mass. 
By combining a calculation of nuclear matrix elements within the framework of the microscopic 
interacting boson model (IBM-2) with an improved calculation of phase space factors, we set limits 
on the average light neutrino mass and on the average inverse heavy neutrino mass (flavor violating 
parameter). 
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The process Oi^/3/3 in which a nucleus X is transformed into a nucleus Y with the emission 
of two electrons and no neutrinos, z^n —^ z+ 2 ^n- 2 +‘^Q~, is of fundamental importance for 
determining the Majorana or Dirac nature of the neutrino and conhrming a non-zero value 
of its mass as established by neutrino oscillation experiments [l-3|, what constitutes physics 
beyond the standard model. The half-life for this process can be written as 
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where Gqi, is a phase space factor (PSF), is the nuclear matrix element (NME), and 
/ contains physics beyond the standard model through the masses m* and elements Uei 
of the mixing matrix of the neutrino (or other hypothetical particle beyond the standard 
model). We have recently (i) introduced a new method j4|, the microscopic interacting boson 
model, IBM-2, to calculate the NME in a consistent way for all nuclei of interest, and (ii) 
improved the calculation of the phase space factors (PSF) by solving the Dirac equation 
for the ongoing electrons in the presence of a charge distribution and including electron 
screening ^ . In this letter, we present results of a calculation that combines the NMEs and 
the PSFs to half-lives. By comparing with current experimental limits we then set limits on 
neutrino masses and their couplings. 

Starting from the weak Lagrangean, £, one can derive the transition operator inducing 
the decay, which, under certain circumstances, can be factorized as T{p) = H{p)f{mi, Uei), 
where p = is the momentum transferred to the leptons 6l-l8|. The transition operator 
H{p) has the form 


Hip) = [-h^ip) + h^^ip) 


O^T) * O^n ' 


+h^ip)SL'] ■ 


( 2 ) 


The factors h^’^'^’'^ip) are given by h^’^'^’'^ip) = r;(p)/i^’'^^’^(p), where v{p) is called the 
neutrino ’’potential” and h{p) are the form factors, listed in Ref. 8j. This form assumes 
the closure approximation which is expected to be good a approximation for 012(3(3 decay 
since the neutrino momentum is of the order of 100 MeV/c while the energy scale 
of the nuclear excitations is 1 MeV, and all multipoles in the intermediate odd-odd nucleus 
contribute to the decay. (Conversely, the approximation is not expected to be good for 
21/(3(3 decay, where the neutrino momentum is of order 2 MeV/c, and only I"*" and O’*" states 
in the intermediate odd-odd nucleus contribute to the decay). The hnite nucleon size is 
taken into account by taking the coupling constants momentum dependent and short range 
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correlations (SRC) are taken into acconnt by convoluting v{p) with the correlation function 
J{p) taken as a Jastrow function. The functions /(m*, Uei) and H{p) depend on the model 
of decay. We consider here explicitly two cases: (i) the emission and reabsorption of a 
light {miight 1 keV) neutrino; (ii) the emission and reabsorption of a heavy {rriheavy ^ 1 
GeV) neutrino. For scenario (i), the function / can be written as 

/ = ■^^, {my) = ^ {Uekf ruk, (3) 

rUe 

k=light 


where U is the neutrino mixing matrix. The average neutrino mass is given in terms of mixing 
angles and phases [n| and is constrained by atmospheric, solar and neutrino oscillation 
experiments. The potential v{p) for this case is v{p) = 27r~^[p{p + where A is the 

so-called closure energy. For scenario (ii) the transition operator can be written as Th{p) = 
Hh{p)fhijni) Uei), where the index h refers to heavy. The function fh can be written as 

1 \ / 1 \ .9 1 


fh = rup 


rrih 


— ) = 
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k=heavy 




( 4 ) 


The neutrino potential is Vhfp) = 27r“^(merup)”^. The function fh is often written as p 
and called the flavor violating parameter. The average inverse heavy neutrino mass has in 
the past been considered as an unconstrained parameter. However, recently, it has been 
suggested that some constraints can be put on this quantity from large hadron collider 
(LHC) physics and lepton flavor violating processes. The effect of heavy neutrinos on neu¬ 
trinoless double-/? decay has been illustrated within the context of a specihc model as a 
function of the mass of the lightest heavy neutrino in the range 1-500 GeV. 

We have calculated the nuclear matrix elements within the framework of the microscopic 
interacting boson model, IBM-2 13|, in all nuclei of interest. Details of the calculation 
are given in Ref. jj] and in a forthcoming long publication ji^ . Matrix elements 
for light neutrino exchange are shown in Table [I] and Fig. [H where they are compared 
with those calculated with other methods, most notably QRPA Q and ISM |l^ with 
the same (or similar) approximations for the SRC. We note both in Table [I] and Fig. [1] 
a close correspondence between the IBM-2 and QRPA calculations, while the ISM results 
are approximately a factor of 2 smaller than IBM-2/QRPA. (The origin of the difference 
is not completely clear. The three models make different approximations and at different 
levels. A recent combined analysis of Oz//?/? and 2u(3(3 decay seems to indicate that the 
main difference is the size of the model space in which the calculations are done. This is 
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TABLE I. Neutrinoless double-/3 decay matrix elements in IBM-2 with Argonne CCM SRC 

and qa = 1.269, in QRPA with Argonne CCM SRC and gA = 1-254, and ISM with UCOM SRC 
and QA = 1.25. 


A 

IBM-2 

QRPAf. 

ISM^ 

48 

2.28 


0.85 

76 

5.98 

5.81 

2.81 

82 

4.84 

5.19 

2.64 

96 

2.89 

1.90 


100 

4.31 

4.75 


no 

4.15 



116 

3.16 

3.54 


124 

3.89 


2.62 

128 

4.97 

4.93 

2.88 

130 

4.47 

4.37 

2.65 

136 

3.67 

2.78 

2.19 

148 

2.36 



150 

2.74 



154 

2.91 



160 

4.17 



198 

2.25 




Ref. SI 
Ref. [16 1 


substantiated by the observation that the behavior with mass number of all three calculations 
is similar and that they can be reconciled by a simple renormalization). Matrix elements 
for heavy neutrino exchange are shown in Table HIl By combining the matrix elements 
with the phase space factors of Ref. j^, we obtain the expected half-lives shown in Table HTTl 
left, and Fig. [2] for light neutrino exchange and Table ITVl left, for heavy neutrino exchange. 
It should be noted that the combination must be done consistently. If the phase space 
factors of Ref. are used, the nuclear matrix elements of Tables [I] and [IT] must be 
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FIG. 1. (Color online) Nuclear matrix elements for Or^/3/3 decay in IBM-2 compared with 

QRPA [1^ and ISM [l6[. 



FIG. 2. (Color online) Expected half-lives for {ruv) = 1 eV, qa = 1.269. The points for ^^®Te and 
^^®Nd decays are not included in this figure. The figure is in semilogarithmic scale. 

multiplied by g\, that is Mq^, = in Eq. (1). 


Using the experimental upper limits from a compilation of Barabash 18|, the IBM-2 


Q. 


matrix elements of Tables [H and mi and the phase space factors of 15||, we estimate current 
limits on the neutrino mass given in Tables IIIIl right, and Table IIVI right, which are the 
main results of this letter. In Table HVl we give limits both on the flavor violating parameter 
rj and on the average heavy neutrino mass, dehned as {m^^)/mp = where 
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TABLE II. Neutrinoless double-/? decay matrix elements in IBM-2 with Argonne CCM SRC 

and QA = 1.269, and in QRPA with Argonne CCM SRC, gA = 1-25 and intermediate size for the 
model space. 


A 

IBM-2 

QRPAf. 

48 

46.3 


76 

107 

233 

82 

84.4 

226 

96 

99.0 


100 

165 

250 

no 

155 


116 

no. 


124 

79.6 


128 

101 


130 

92.0 

234 

136 

72.8 


148 

103 


150 

116 


154 

113 


160 

155 


198 

104 




Mw = 80.41 ± 0.10 GeV and Mwr is assumed to be Mwr = 3.5 TeV. While the former is 
model independent, the latter depends on the model of left-right mixing jl^ . 

These results are obtained using the free value of the axial vector coupling constant as 
obtained from neutron decay, qa = 1.269. It is known from single {3 decay and 2z//?/? decay 
that Qa is renormalized in nuclei. There are two reasons for the renormalization: (i) the 
limited model space within which the calculation of the NME is done; (ii) the omission of 
non-nucleonic degrees of freedom Since the coupling constant qa appears to the 

fourth power in the life-time, the renormalization effect is non negligible and it will amount 
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TABLE III. Left; Calculated half-lives in IBM-2 for neutrinoless double-/? decay for {m^J) = 1 eV 


and gA = 1-269. Right: 
compilation of Barabash 
from IGEX 


Upper limit on neutrino mass from current experimental limit from a 
181. The value reported by Klapdor-Kleingrothaus et al. [l^, the limit 


20], and the recent limits from KamLAND-Zen 2]| and EXO 


22l| are also included. 


Decay 

T^J^ilO^^yr) 

^l/2,exp(y^) 

{m^){eV) 

^SCa^^^Ti 

0.782 

> 5.8 X 10^2 

< 3.7 

7®Ge^™Se 

1.19 

> 1.9 X 10^5 

< 0.25 



1.2 X 1025f_ 

0.32 



> 1.6 X lO^*^^ 

< 0.27 


0.423 

> 3.6 X 1023 

< 1.1 

96Zr^96MQ 

0.588 

> 9.2 X 102i 

< 8.0 

ioomo^iooRu 

0.340 

> 1.1 X 1024 

< 0.56 

iiopd^iiocd 

1.22 




0.602 

> 1.7 X 1023 

< 1.9 

I24sn^i24xe 

0.737 



i28Te^i28xe 

6.94 

> 1.5 X 1024 

< 2.2 

isoTe^isoxe 

0.355 

> 2.8 X 1024 

< 0.36 


0.512 

> 5.7 X 1024^ 

< 0.30 



> 1.6 X 102'^^ 

< 0.18 

148Nd^l48Sm 

1.79 



isoxd^isogj^ 

0.213 

> 1.8 X 1022 

< 3.4 


3.94 



i®°Gd^^6°Dy 

0.606 



198pt^l98jjg 

2.64 




^ Ref. 
Ref. 
Ref. 
Ref. 


20 


n 


22 | 


to a multiplication of the limits in Table IIIII and IIVI by a factor of 2-4. Details of the 
renormalization procedure, as well as of the calculation of the renormalized matrix elements 
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TABLE IV. Left: Calculated half-lives for neutrinoless double /3 decay with exchange of heavy 


neutrinos for r] = 2.75 x 10 ' and qa = 1.269. Right: Upper limits of \r]\ and lower limits of heavy 


neutrino mass from current experimental limit from a compilation of Barabash 

Q, 


reported by Klapdor-Kleingrothaus et al. 


from KamLAND-Zen 


the limit from IGEX 


211] and EXO 


Q- 


The value 


20(|, and the recent limits 


are also included. 


Decay 



|r/|(10-7) 

(m^J(GeV) 

^®Ca^^®Ti 

0.096 

> 5.8 X 10^2 

< 3.54 

> 0.73 

76Ge^76se 

0.190 

> 1.9 X 1025 

< 0.275 

> 9.4 



1.2 X 1025_^ 

0.346 

7.5 



> 1.6 X 1025_^ 

< 0.300 

> 8.6 


0.070 

> 3.6 X 1025 

< 1.22 

> 2.1 

96Zr^96MQ 

0.025 

> 9.2 X 1021 

< 4.56 

> 0.6 

IOOMo^IOORu 

0.012 

> 1.1 X lO^'^ 

< 0.285 

> 9.1 

llOpq^llOcd 

0.044 





0.025 

> 1.7 X 1025 

< 1.06 

> 2.5 

I24sn^i24xe 

0.089 




i28Te^i28xe 

0.846 

> 1.5 X 1024 

< 2.07 

> 1.2 

isoTe^isoxe 

0.042 

> 2.8 X 102“^ 

< 3.38 

> 7.6 


0.066 

> 5.7 X 1024_^ 

< 0.296 

> 8.7 



> 1.6 X 1025_f_ 

< 0.177 

> 14.6 

148Nd^l48Sm 

0.048 




isoxd^isogj^ 

0.006 

> 1.8 X 1022 

< 1.58 

> 1.6 


0.132 





0.022 




198pt^l98jjg 

0.063 





Ref. 
Ref. 
" Ref. 
Ref. 


n 


m 

21 


22 ] 


NME, will be given in a forthcoming longer publication |l^. The question of whether or 
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FIG. 3. (Color online) Current limits to {my) from CUORICINO 


KamLAND-Zen 


2M, IGEX 


20], NEMO-3 


21(1, and EXO 2^ and IBM-2 nuclear matrix elements. The value of Ref. 


is shown by X. It is consistent only with nearly degenerate neutrino masses. The figure is in 
logarithmic scale. 


not 01/13(5 matrix elements should be renormalized as much as 2u(3(3 matrix elements is the 
subject of much debate. In 2i/(3(3 only 1+ and O’*" states in the intermediate odd-odd nucleus 
contribute to the decay, while in 0i/(3(3 all multipoles play a role. In this letter we do not 
dwell on this question, but rather present results with the unrenormalized value qa = 1.269, 
summarized in Fig. |3l From this hgure, one can see that in the immediate future only 
the degenerate region can be tested by experiments and that the exploration of the inverted 
region must await much larger (> Iton) experiments, especially if qa in 0i/(3(3 is renormalized 
as much as in 2i/(3/3 decay. From the same hgure, one can also see that even the one-ton 
experiments will not be able to reach into the normal hierarchy. 

This work was performed in part under the USDOE Grant DE-FG02-91ER-40608 and 
Fondecyt Grant No. 1120462. 
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